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ABSTRACT 

The heating of the lower solar corona is examined using numerical simula- 
tions and theoretical models of magnetohydrodynamic turbulence in open mag- 
netic regions. A turbulent energy cascade to small length scales perpendicular to 
the mean magnetic field can be sustained by driving with low-frequency Alfven 
waves reflected from mean density and magnetic field gradients. This mechanism 
deposits energy efficiently in the lower corona, and we show that the spatial dis- 
tribution of the heating is determined by the mean density through the Alfven 
speed profile. This provides a robust heating mechanism that can explain ob- 
served high coronal temperatures and accounts for the significant heating (per 
unit volume) distribution below two solar radius needed in models of the origin 
of the solar wind. The obtained heating per unit mass on the other hand is 
much more extended indicating that the heating on a per particle basis persists 
throughout all the lower coronal region considered here. 

Subject headings: Sun: corona — MHD — turbulence 



1. Introduction 

The mechanism by which the solar corona is heated is a fundamental problem in astro- 
physics and plasma physics, and one that remains incompletely explained. An acceptable 
mechanism transports energy to the corona, heating it to temperatures greatly in excess 
of chromospheric values. This mechanism accelerates the solar wind and establishes the 
boundary conditions for the entire plasma heliosphere. There are several constraints on an 
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acceptable coronal heating model. Spectrometers (Kohl et al. 1995; Kohl et al. 1997) show 
that proton temperatures in polar coronal holes rise to several 10 6 K at heliocentric distances 
of 1.5 to 2 solar radii (R s ) and the bulk wind speed reaches 200km/s. Remote sensing ob- 
servations (Grail et al. 1996) indicate wind speeds up to 600 km/sec at 4R S . Fluid models 
of the solar wind (Habbal et al. 1995; McKenzie et al. 1995; Evje & Leer 1998) support the 
conclusion that wind acceleration at such low altitudes requires heating low in the corona. In 
the fluid models an ad hoc heating function Q(r) is assumed, which represents the energy de- 
position per unit volume as Q = Q exp [— ^^p] , where r is heliocentric distance, r = R s 
7 x 10 5 km, L is a prescribed dissipation lengthscale, and Q is chosen so that the energy flux 
has the required value of about 5 x 10 5 erg cm" 2 s _1 (Withbroe & Noyes 1977). Despite the 
fact that there is no direct observational evidence for the exponential heat function, it has 
been long known to provide (Holzer & Axford 1970; Kopp & Orrall 1976; Hammer 1982) 
rapid acceleration of the wind, and with L m 0.2-0.5 R s it accounts for many observed fea- 
tures of the corona (McKenzie et al. 1995; Habbal et al. 1995; Axford & McKenzie 1997). 
Although some possible sources of heating are proposed, in those models there is no physical 
mechanism for explaining this heating function and, being an ad hoc quantity, the exact ex- 
plicit exponential function expression is perhaps irrelevant and the main point to be remarked 
is the requirement of strong heating per unit volume in the first solar radius of the corona. 
Yet other two-fluid (electrons and protons) models (Tu & Marsch 1997; Marsch & Tu 1997) 
have included some mechanisms for the heating, based on the requirement of high fre- 
quency waves launched from the Sun, but those high frequencies (of the order of KHz) 
waves have yet remained unobserved in the coronal base where the energy is injected. In 
another previous approach, Hollweg 1986; Hollweg & Johnson 1988; Isenberg 1990 conjec- 
tured that waves are launched from the Sun at low frequencies but magnetohydrodynamic 
(MHD) turbulent-like cascades are assumed to transfer the energy to higher frequencies, 
where ion cyclotron mechanisms can be invoked to dissipate it. This idea, recently em- 
ployed by Li et al. 1999 and reviewed by Cranmer 2000 involve phenomenological turbu- 
lent rates (being called of the Kolmogorov and Kraichnan type) but two important prob- 
lems are faced. First, a cascade in wave frequencies is assumed, which implies a turbu- 
lent cascade proceeding in the parallel (to the mean field) wavenumber direction and a 
linear dispersion relation to persist between parallel wavenumber and frequency. MHD 
simulation studies show on the other hand that in the presence of a strong mean mag- 
netic field, like in the corona, the parallel wavenumber cascade is almost completely su- 
pressed (Shebalin et al. 1983; Oughton et al. 1994; Matthaeus et al. 1998) in favor of the 
perpendicular cascade. This is precisely the kind of anisotropy which is well documented 
(Armstrong et al. 1990) for the lower coronal density. 

Second, the necessity of counterpropagating fluctuations to sustain the turbulence 
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(Kraichnan 1965; Dobrowolny et al. 1980) is not considered in the models and a WKB ap- 
proach is implicitly assumed. In those models, the turbulent cascade was not the central 
focus of study, but we believe that a useful perspective can be gained into the heating prob- 
lem by further study of this issue. In this present work, instead of assuming a parallel wave 
cascade, we describe ab initio the conditions under which a strong (perpendicular) turbu- 
lent cascade can actually be developed and sustained by low frequency driving at the base. 
Reflections from density and mean field gradients provide the necessary counterpropagating 
fluctuations to sustain the turbulence and this effect is fully taken into account in our model. 
We will employ both direct numerical simulation and turbulent phenomenologies to obtain 
the heating rate distribution. 

We show that the heating distribution profile is related to the coronal background 
magnetic field and density profile, i.e., the heating and the density profile are not independent 
quantities. On the one hand, we find that the heating per unit volume is significant mostly 
within the first solar radius, with a dissipation length scale (defined as the length in which the 
heating decays by an order of magnitude) of a fraction of solar radius. On the other hand the 
heating per unit mass varies much less dramatically and is significant through the complete 
domain. For a complete model, appropriate kinetic processes will need to be identified to 
absorb the small perpendicular scale cascaded energy. Candidates for small scale damping are 
plentiful, and include oblique kinetic Alfven waves, collisionless reconnection and associated 
heating at current sheets and nonlinear beam instabilities (Leamon et al. 1998). The kinetic 
physics of heating is beyond the scope of the present paper. The main purpose of this work 
is to show that a heating distribution emerges naturally from a consistent MHD turbulent 
cascade driven by low frequency fluctuations at the coronal base. This result may bring us 
closer to understand the puzzle of coronal heating and the origin of the solar wind. 



2. Model 

Most coronal heating models derive their energy supply from either wave power or qua- 
sistatic field line motions, originating in the photosphere or chromosphere. These differ 
mainly in their characteristic timescales. A general scenario for wave production in open 
magnetic regions has been proposed (McKenzie et al. 1995) in which activity in the net- 
work regions, driven by convective motions, provides a source of upward traveling Alfven 
waves. Higher frequency waves are easily transported upward, whereas low-frequency 
waves experience greater reflection and transmit less readily into the corona (Parker 1965; 
Heinemann & Olbert 1980; An et al. 1990; Zhou & Matthaeus 1990; Velli 1993). However 
even in the limit of sharp gradients or discontinuities, there are indications (Hollweg 1981) 
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that a substantial energy flux can reach the coronal base (Kudoh & Shibata 1999). 

Here we examine the following scenario (Matthaeus et al. 1999): MHD-scale ( low- 
frequency ) fluctuations are excited below the base of the corona, on transverse length 
scales that are characteristic of the chromospheric network, and propagate upwards along 
the mean magnetic field B. A portion of this wave flux enters the corona. Reflections from 
large-scale inhomogeneities within the lower corona produce counter-propagating fluctua- 
tions (Heinemann & Olbert 1980; Zhou & Matthaeus 1990) a situation that permits strong 
nonlinear MHD couplings (Dobrowolny et al. 1980). The couplings drive a perpendicular 
(to B) turbulent cascade, producing small scale reconnection events that couple to kinetic 
processes at small scales and oblique wavevectors (Leamon et al. 1998). 

In our model, we simulate the propagation of low-frequency MHD fluctuations in a 
medium whose properties are chosen to represent an open magnetic field region (like a polar 
coronal hole) in the lower corona. We consider a volume section which extends in the radial 
direction a distance L s of order 2R S above the surface. A sketch of such a region is shown 
in Fig. 1. 

Our goal is to compute the dynamics of transverse fluctuations of magnetic field b and 
velocity v as influenced by specified large-scale inhomogeneities (mean magnetic field and 
density gradients) as well as by local nonlinear couplings. To do this, in a medium where 
there is a strong mean magnetic field, we use the reduced MHD (RMHD) approximation 
(Strauss 1976; Montgomery 1982; Zank & Matthaeus 1992). This approach is appropriate 
for fluctuations in a nearly incompressible plasma, where the plasma beta (ratio of thermal 
to magnetic pressure) is low (Zank & Matthaeus 1992). We employ the Elsasser variables 
z± = v ± b (where the magnetic field is written in units of b/ y/Anp, with p=density) to 
represent downward (+) and upward (-) propagating fluctuations. Fluctuations depend on 
the radial coordinate r, the transverse spatial coordinates (including an areal expansion 
factor A(r)) and on time t. Mean quantities like density p and magnetic field B are assumed 
to depend only on r. Consistency with the RMHD assumptions (Zank & Matthaeus 1992) 
requires d r <C Vj_, Va 3> z±, and Vj_ • z± = 0, where Va = B / y/Airp is the Alfven speed 
which gives the propagation velocity of fluctuations in the linear limit. 

In the presence of a non-uniform Va, the nonlinear RMHD equations can be written 
(Zhou & Matthaeus 1990) as 

^r + V A ^ = -Riz+ + R 2 z- - z + • V ± z_ + rzViz_ 
Ot Or 

-^--Va-^- = Riz- - i? 2 z+ - z_ • V ± z+ + V Vlz + (1) 
ot Or 

where r\ is the resistivity (assumed equal to the viscosity), and Ri(r),R 2 (r) are specified 
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reflection rates due to the inhomogeneous Va(t). For a symmetric radial geometry (A(r) = 
r 2 ), B = B (r /r) 2 (with B the magnetic field at r = R s ), the reflection rates are R\(r) = 
and i?2( r ) — + v - - m the case of a general geometry with cross section area A(r), 
where A(r) = A B /B(r), the reflection rates are R^r) = \^ and R 2 (r) = + 
The previous expressions are valid for a section with transverse dimensions small compared 
to the longitudinal size and close to the polar region in the Sun (dependence of the mean field 
B with the latitudinal angle should be assumed to maintain the free divergence condition, 
V-B = 0). 

A Chebyshev-Fourier representation is chosen to solve the equations (1) using a di- 
rect numerical simulation of the pseudospectral type. Transverse coordinates are periodic, 
while non periodic boundary conditions are imposed on the bottom and top sections of the 
volume. In order to model the input of velocity fluctuations at the base of the coronal 
domain, we refer to observations of non-thermal velocity amplitudes in the upper transi- 
tion region of the Sun. These amplitudes are reported to range from 20 to 55 km/sec 
(Chae et al. 1998; Doyle et al. 1998; Hassler et al. 1990). The system of nonlinear equa- 
tions (1) is driven by imposing a time dependent oscillatory pattern of velocity and magnetic 
field at the bottom surface, with a transverse lengthscale corresponding to supergranule diam- 
eters of 13-32 Mm at the coronal base (Hagenaar et al. 1997) and the typical inter-network 
spacing of 30 Mm (Axford & McKenzie 1997). Transverse dimensions of the simulation box 
are chosen to include about 4 supergranules. In the bottom panel of Fig. 2 the boundary 
conditions imposed on the form of waves at the bottom of the domain are shown. The 
top panel of Fig. 2 shows a cross section of the volume above in the coronal simulation 
(corresponding to one of the simulations to be described below). Numerical resolution of 
128 x 128 x 33 Fourier +Chebyshev grid points are employed in this simulation, for a macro- 
scopic Reynolds number of order 600. The fluctuating fields are highly structured on the 
transverse directions as shown by concentrations of current density in the form of sheet-like 
small scale structures. 

Driving is at a low frequency of / = 0.6 /t a, where tA = Rs/Va is a vertical Alfven 
crossing time based on the Alfven velocity at the base and a distance of a solar radius. This 
correspond to driving wave periods of the order of 1000 sec. 

Open boundary conditions are assumed at the top surface. If reflections were absent, 
this would lead to a wavetrain of upward propagating fluctuations that would escape contin- 
uously from the upper boundary without depositing any energy in the system. The system 
would dynamically relax toward a maximal cross helicity state (Dobrowolny et al. 1980; 
Grappin et al. 1983; Ting et al. 1986) in which downward propagating fluctuations vanish 
and there would be no incompressible MHD nonlinearities and hence no cascade. However, 
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the first two terms on the right hand sides of Eqs. (1) are associated with gradients of the 
wave speed, and these lead to partial reflection of the excited waves. This results in a state 
with both upward and downward type fluctuations and thus, sustained turbulence. We have 
presented in a previous numerical study (Dmitruk et al. 2001a) the conditions under which 
an efficient MHD turbulent perpendicular cascade can be sustained; what is required is oc- 
currence of reflections and non-propagating structures (controlled by the type of boundary 
conditions imposed). The efficiency of this mechanism on dissipating the injected energy has 
been demonstrated by both phenomenologies (Matthaeus et al. 1999; Dmitruk et al. 2001b) 
as well as direct numerical simulations (Dmitruk et al. 2001a; Oughton et al. 2001). Effi- 
ciencies (ratio of turbulent dissipation rate by the energy injection rate) averaging between 
10 and 40 % can be obtained. Even very weak reflections are capable of sustaining the 
turbulence and provide effective dissipation. 

In the next section we obtain the heating distribution profiles resulting from this tur- 
bulent state. 



3. Heating profiles 

We performed numerical simulations employing several different mean magnetic field 
and density radial profiles. The values at the bottom boundary (Axford & McKenzie 1997; 
Feldman et al. 1997) have been fixed to B = 9 Gauss, and a number density of n = 
3.2 x 10 8 cm -3 , which gives an Alfven speed at the base of 1100 km/s. With the assumed 
amplitude of fluctuations of the order of 20 — 55 km/s, this gives an average energy input flux 
per unit area, F^ ~ PoVa < z 2 _ — z\ > of order 5 x 10 5 erg cm" 2 s _1 corresponding to the 
Withbroe & Noyes 1977 value. These values have to be taken as order of magnitude quan- 
tities but a different choice would not invalidate the general argument we want to present 
here. In our first simulation we consider a symmetric radial mean magnetic field profile 
B{r) = B (R s /r) 2 . The mean density follows from the hydrostatic equilibrium approxima- 
tion (An et al. 1990; Velli 1993; Nakariakov et al. 2000): dp/dr = -rriing, with p = thermal 
pressure, n = number density, rrii the proton mass and g = gravity. In an isothermal corona, 
this equation takes the form dn/dr = —oct nR s /r 2 where ax = pGM / R g TR s (// mean molec- 
ular weight, G gravitational constant, M solar mass, R g gas constant and T temperature). 
The quantity a r 12 for the Sun and a fully ionized hydrogen coronal plasma at T = 10 6 
K. The number density distribution is then n(r) = hq exp[a<r(i? s /r — 1)] and the density is 
p{r) = rriin(r). We consider two such density profiles, shown in the top left panel of Fig. 3 
with continuous and dotted lines, corresponding to values of ax = 12 and «y = 6 (represent- 
ing a hotter corona). This gives two different density profiles, with the hotter corona having 
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the shallower profile. The corresponding Alfven speed profiles Va(t) = B(r)/^Anp(r) are 



shown on the top right panel of Fig. 3. When the system reaches a statistically steady state 
the turbulent dissipation rate e (equal to the cascade rate) is computed from the fluctuating 
fields z_, z + from which the heating per unit volume is obtained as Q = pe. 

This quantity is plotted as a function of r in the bottom left panel of Fig. 3. It is 
apparent that the dissipation per unit volume is peaked near the model coronal base. It is 
important to note that Q(r) decreases faster for the model with the steeper density profile. 
The radial distance Ar at which Q decreases by an order of magnitude with respect to the 
coronal base value Qo is 0.36i? s and 1.1 R s for the two profiles in Fig. 3. This would be 
the analog of the dissipation length in the ad-hoc exponential heating distribution profiles. 
The bottom right panel on Fig. 3 shows the heating per unit mass, a useful quantity to 
consider the heating in a per particle basis. Over our entire domain this quantity maintains 
a relatively constant value. Extended heating per unit mass fits the apparent requirement 
of a hot (non adiabatic) solar wind and has also emerged in empirical models of radiation 
loss and other heating models (Mullan and Cheng 1994; Gibson 1973). 

A third simulation has been performed employing a super radial geometry, where the 
mean magnetic field profile has been taken as 



Profiles of this type have been considered in Hollweg 2000 and similar super expansion 
profiles have been used by McKenzie et al. 1995; Evje & Leer 1998; Li et al. 1999 in their 
solar wind models. We have taken f max = 6 which gives a boundary value of Bo = 9 Gauss. 
An observationally based number density is assumed in this case (Feldman et al. 1997), 



The density profile and resulting Alfven speed profile are shown in the top panels of Fig. 3 
(dashed line). The obtained heating per unit volume and heating per unit mass are shown 
in the bottom panels of Fig. 3 (dashed line). The result confirms the tendency shown for the 
earlier cases: this profile, with the fastest decreasing density, also has the fastest decreasing 
heating per unit volume. The decay length for this case is Ar = 0.2R S . The heating per 
unit mass remains at a relatively constant value. 

The relevant results to be extracted from these simulations for different Alfven speed 
profiles is that the heating per unit volume has significant values only within the first solar 
radius of the lower corona and that the radial profile of the heating is controlled by the radial 




Gauss 





(3) 



- 8- 



profile of the density. As it will be more clear from the phenomenological analysis that follows 
in the next section, it is actually through the Alfven speed profile and its gradient that this 
behavior is manifested. 

Whether the radial profile of the heating per unit volume is exactly exponential or not 
is, we believe, irrelevant and we conjecture that a steep heating profile, like the ones we 
obtained here, would work as well to explain the fast acceleration of the solar wind. 

Fig. 3 contains the essence of the present results, that a reasonable and efficient heating 
can be obtained from the MHD perpendicular turbulent cascade driven by low frequency 
waves injected from the bottom. The properties of this heating are controlled by both the 
amount of energy injected into the system and the background profiles of magnetic field 
and mean density. To further illustrate this idea and, in some limiting case, to obtain an 
analytical expression for the heating profile in the low region of the volume considered here, 
we present in the following section a phenomenological approach to the previous model Eqs 
(!)• 

4. Heating profile from a phenomenological model 

A nonlinear phenomenology (Dmitruk et al. 2001b) simplifies the study of the wave- 
driven RMHD heating model. We modify Eqs. (1) as follows: the independent variables, 
now designated as Z±(r), are treated as one dimensional wave amplitudes with the same 
linear transport terms (LHS and reflection terms) as in Eqs. (1), but with the nonlin- 
ear and dissipative terms on the RHS replaced by the nonlinear models Z T \Z±\/2\j_(r). 
This model entirely suppresses the transverse structure, with the strength of the nonlin- 
ear cascade effects represented through a single correlation lengthscale Aj_(r). Previous 
investigations (Dmitruk et al. 2001b) have shown that this phenomenology portrays many 
of the same physical effects as wave-driven models based upon the full RMHD equations 
(Dmitruk et al. 2001a; Oughton et al. 2001). We expect this model to correspond to high 
Reynolds number fully developed RMHD turbulence since the amplitudes Z_ and Z + are 
strongly coupled in the adopted forms for the nonlinear terms. We extracted a heating 
function Q(r) from a numerical solution to this phenomenological model, using the same 
parameters as those employed in the full nonlinear RMHD simulations and assuming a cor- 
relation lengthscale X±(r) linearly increasing with radius, with a value at the coronal base 
set to correspond to the inter-network spacing of 30 Mm. The derived Q(r) is illustrated in 
Fig. 4 as the thin line. The direct numerical solution of the previous section is shown as 
the thick line. This corresponds to the simpler hydrostatic density and 1/r 2 mean magnetic 
field profile (similar results are obtained for the super expansion geometry). Despite some 
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differences on the far region, where dVA/dr approaches and the heating is very low, the 
agreement is pretty good, considering the extreme simplification involved in the phenomeno- 
logical treatment. Of course, the phenomenology would not be useful for studies of the 
turbulent cascades and spectra, but it gives the right order of magnitude for the averaged 
quantities such as the dissipation (heating), when this quantity is not extremely low. The 
interesting issue of how the correlation length and the strength of the cascade depend on the 
external parameters warrants further study comparing the phenomenology with the full non 
linear direct numerical simulation. 

In Fig. 5 we explore the phenomenological heat function Q obtained by varying the 
correlation lengthscale at the coronal base, with Aj_ = 1 corresponding to 30 Mm. The limit 
Aj_ — > corresponds to stronger local turbulence since the implied eddy turnover time \±/z_ 
becomes much smaller than the wave transit time R s /Va- Fig. 5 shows that Q increases as 
Aj_ — > 0, while maintaining the fast radially decreasing behavior. 

In this limit, an analytical asymptotic expression can be obtained. Although limited 
in its validity, this analysis serves to illustrate better the properties of the heating in this 
model. We begin by formulating a flux balance equation, based upon Eq.(l). This takes the 
form dF/dr = —A(r)Q, where the net upward energy flux F = pA Va(Z 2 _ — Z\) and A is 
the cross section area. Note that in geometrical optics (or "WKB" theory for noninteracting 
waves) (Jacques 1977) the flux balance formally requires that Q = 0. When fluctuations are 
damped, Q ^ (Hollweg 1986). Here the crucial relation that connects the heating to the 
Alfven speed profile emerges by identifying Q with a physically correct turbulent dissipation 
rate per unit volume at distance r, namely Q = p e, where the dissipation rate per unit mass 
e is expressed in the phenomenological model as e ps (Z\ \ZJ\ + Z 2 _\Z+\)I \±. We now proceed 
to form an asymptotic strong turbulence limit as follows. Letting A^ — > 0, one can conclude 
(Dmitruk et al. 2001b) that the downward flux is small, \Z + \ <C |Z_|, and that in a steady 
state this limit is characterized by |Z + |/Aj_ = \dVA/dr\. Thus e ps Z 2 _ |Z + |/Aj_ and we obtain 
the asymptotic result that the dissipation per unit length can be related to the Alfven speed 
profile and the energy flux: A Q ps pA Z 2 __\dV aI dr\ ps V^^Va/ dr\F . Thus, the asymptotic 
flux balance equation can now be written as dF/dr ps — V^^Va/ dr\F , whose solution is 

Vao/Va if r < r m 

F Q (VaJVaJ 2 V a /V Ao if r>r m U 

where F Q is the net flux at the base and r m is the radial distance at which Va reaches its 
maximum value VA m . 

This allows us finally to obtain an explicit expression for the heating deposition 

^ f F Ao {Ml A) \dV A /dr\ (VaJVD if r < r m 
W ~ I F Aq (A /A) \dV A /dr\ {VaJVIJ if r > r m ^ 



-10- 



where Fa = F /A is the energy flux per unit area at the base (i.e, the Withbroe & Noyes 1977 
value) . 

For the particular case of the A = r 2 geometry and the isothermal atmosphere used 
before, with r m 3R s , this relation predicts the top dot dashed curve of Fig. 5. As seen in 
this figure, the phenomenological solutions for different Aj_ values approach this asymptotic 
case as Aj_ tends to 0. Already for A^ = 1/10 the asymptotic solution would be a good 
approximation. This corresponds to correlation length scales of order 3000 km. Although 
much less than the typical inter network distances, observations do not rule out correlation 
scales of this magnitude. Eq. (5) also shows that in this limiting case the value of the heating 
goes to at points where dVA/dr = 0, a limitation which is only apparent at corresponding 
very low values of Q in the numerical solution (a similar result applies for the super expansion 
profile presented before, where the point at which dVA/dr = is r m ps 1.6R s ). 

For this asymptotic limit, eq. (5) explicitly shows that the spatial distribution of tur- 
bulent heating is directly connected to the radial profile of the Alfven speed and its absolute 
value depends on the energy input flux. This illustrates the connection already pointed out 
through the direct numerical simulations results. 

5. Conclusions 

We have presented a model for MHD turbulence driven by low-frequency waves im- 
posed at the lower boundary. We applied this model to an open region in the lower corona, 
adopting reasonable profiles for magnetic field and density. Both direct numerical simula- 
tions and phenomenological studies have been performed and compared. Results show that 
turbulence can be sustained, due to reflection of upward traveling waves by Alfven speed 
gradients, producing a perpendicular turbulent cascade and efficient heating. The energy 
dissipation per unit volume is significant mostly in the lower part of the model corona with 
a dissipation length scale of a fraction of solar radius. The obtained heating profile Q is sim- 
ilar to profiles commonly used in solar wind acceleration models. However the more general 
result emerges, for strong turbulence, that the heat function Q(r) should be related to the 
Alfven speed profile Va{t). This shows that the heating distribution is determined by the 
large-scale magnetic field and density profiles. In this way one sees that the scale-height of 
the heating per unit volume is directly related to the density scale-height, thus explaining 
the spatial confinement of Q as a consequence of the rapid density decrease in the lower 
corona. However, the deposition of heat per unit mass Q/ pis much more extended, perhaps 
well beyond the lower coronal region considered in the present model. These conclusions 
do not rule out other models, such as direct cyclotron absorption of high-frequency waves 
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(Axford & McKenzie 1997; Tu & Marsch 1997) but they imply that it is possible to pro- 
vide the proper kind of heating per unit volume using a low-frequency wave-driven MHD 
turbulence mechanism. 

Research supported by NASA (NAG5-7164, NAG5-8134 SEC Theory Program), NSF 
(ATM-9713595, ATM-0296113), and UK PPARC (PPA/G/S/1999/00059). 
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Fig. 1 - Sketch of the coronal region considered in this model 

Fig. 2 - In the bottom panel the wave forcing applied at the bottom boundary is illustrated 
with the imposed transverse fluctuating field z_ (arrows) overlaying the current density 
(grayscale). In the top panel a cross section at a fraction 0.2R S above the base is shown 
with transverse fluctuating field (arrows) over the current density (grayscale) obtained in 
the numerical simulation. Note the small scale structure formation. 

Fig. 3 - The panel at the top left shows three different number density profiles assumed for 
the simulations. The continuous and dotted line correspond to the hydrostatic model for 
different values of the density scale height. The dashed line correspond to an observationally 
based density (see text). The top right panel is the corresponding Alfven speed profile, with 
1/r 2 geometry for the hydrostatic density and super expansion geometry for the observational 
density. The bottom left panel shows the result of the heating per unit volume from the 
simulations. The bottom right panel is the heating per unit mass. 

Fig. 4 - Comparison between the heating per unit volume obtained with the direct numerical 
simulation of Eqs.(l) (thick line) and the numerical solution of the phenomenological model 
(thin line), for the case of an hydrostatic density profile and 1/r 2 geometry. 

Fig. 5 - Different solutions for the heating per unit volume of the phenomenological model 
(thin lines) for values of the correlation length A^ = 1,0.5,0.1,0.02 (in units of the inter 
network length) which approach the asymptotic solution when A^ — > (Eq.(5) in the text) 
shown as the dot dashed line. 
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